Does the 2D Hubbard Model Really Show d-Wave Superconductivity? 
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In a recent Letter |l[ , Maier et al concluded that they 
found, based on the dynamical cluster as well as non- 
crossing approximations on Green functions, a stable 
pure d-wave solution with off-diagonal long-range order 
(ODLRO) in the intermediate to strong coupling regime 
for the doped repulsive 2D Hubbard model, and given 
an estimation of maximum Tc being \bQK for doping 
X « 0.2. 

This is a surprising yet too strong statement. The au- 
thors in Ref. |^ might be unaware of a theorem, proved 
by Su and Suzuki ||^ two years ago, that the 2D Hub- 
bard model with narrow bands (including next nearest- 
neighbor hopping, etc.) does not exhibit d-wave pairing 
long-range order (LRO) at any nonzero temperature for 
both repulsive and attractive Coulomb interactions and 
for any electron fillings. This rigorous result has been 
generalized to cover other models (see, e.g. [^). There are 
also a few rigorous results showing the absence of s-wave 
or other-wave pairing LROs for the 2D Hubbard model 
at finite temperatures ^J^. It seems that the statement 
in Ref. [Q contradicts these mathematically rigorous re- 
sults. On the aspect of numerical simulations, although 
some quantum Monte Carlo (QMC) calculations show 
tendencies favoring d-wave superconductivity , no def- 
inite sign of LRO has been detected, as pointed out by 
many authors (see, e.g. 0-^). Even though the long- 
tailed enhancements in d-wave pairing correlations near 
half-filling or d-wave-like pairing fluctuations at low tem- 
peratures were observed in QMC studies, none con- 
cluded from QMC results that the d-wave superconduc- 
tivity with ODLRO is clearly obtained. In fact, the QMC 
results, owing to the well-known limitations, though still 
controversial, are not incompatible with the rigorous re- 
sults. There indeed exist analytic and numerical works 
(e.g. QMC) pO|-p^, in completely opposite to the state- 
ment claimed in Ref. |^], showing that the 2D Hubbard 
model (with a next nearest-neighbor hopping integral) 
does not exhibit any definite sign of s-wave and d-wave 
superconductivity. Why are the consequences gained in 
Ref. [Q inconsistent with the rigorous results? 

Now, let us look at what was really observed in Ref. 
101 . First, the method used in Ref. is only approxi- 
mate. They used the coarse grained Green function to 
mimick the true Green function. According to the defini- 
tion, as the coarse grid number Nc tends to infinity, the 
coarse grained Green function recovers the true Green 
function, whereas in the realistic presentations they use 
just the result oiNc = 4 to draw final conclusions without 
giving a convincing comparison for different N^s, apart 



from just stating the systematic errors on the order of 
Second, one may notice that the statement of 
the appearance of a pure d-wave superconductivity with 
ODLRO comes from such an observation that the anoma- 
lous coarse grained Green function vanishes at the points 
(0,0) and (tt, tt) but is finite at (0, tt) and (7r,0). De- 
spite how realiable the approximate method is, one can- 
not, just from this simple observation, conclude that a 
pure d-wave superconductivity is uncovered. For exam- 
ple, a mixing state like d^2_y2 + idxy cannot be ruled 
out. Third, the transition observed in Ref. j^] was specu- 
lated as a possibility of a Kosterlitz-Thoulcss (KT) phase 
transition. As is well-known, the correlation function in 
the KT phase decays in a power law, while the existence 
of ODLRO implies, as discussed by Yang that the 
correlation function becomes a nonzero constant as the 
spatial separation between the pairing operators tends 
to infinity. In their seminar article pH , Kosterlitz and 
Thouless also pointed clearly out that this type of phase 
transition cannot occur in a superconductor. As a re- 
sult, one cannot again just from the approximate results 
conclude the occurrence of ODLRO in the 2D Hubbard 
model at low temperatures. 

In short, the results presented in Ref. are not suffi- 
cient to allow for drawing the conclusion like the authors 
did, at least in a conventional sense. As a matter of 
fact, it appears that the transition they observed might 
be a kind of topological phase transition with a power 
law decay (i.e. possessing a quasi-LRO), not a genuine 
superconducting transition because the latter possesses 
a true LRO. This reinterpretation makes the observation 
in Ref. consistent with the existing rigorous results, 
namely, the 2D Huabbrd model does not show d-wave 
pairing LRO at finite temperatures. 
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